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ABSTRACT 



The time-dependent mixed-layer depth and temperature 
response of a two-layer ocean model to a moving hurricane 
model was investigated. The hurricane model was that of 
Elsberry, Pearson, and Corgnati (1974). The maximum hurri- 
cane wind was bounded at 50 m/sec in order to isolate ocean 
reactions from hurricane variance. The ocean model was 
based on that of Fraim (1973) which was derived from the 
earlier mixed-layer models of Kraus and Turner (1967) and 
Denman (1973). The ocean model was an idealized, quasi- 
three-dimensional model designed to examine the effects 
along the path of movement due to the passage of a moving 
hurricane. Results indicated that slower-moving storms 
engendered the greatest change in mixed-layer depth and 
temperature. Away from the center a combination of entrain- 
ment mixing and convection was responsible for the majority 
of the depth and temperature changes. Near the eye, in the 
upwelling zone, a combination of entrainment mixing, convec- 
tion, and advection caused temperature change while advec- 
tion alone dominated layer depth change. 
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I. INTRODUCTION 



There exists a large body of evidence which demonstrates 
that there is a complex interaction between a hurricane and 
the ocean. For example it appears that sea-surface tempera- 
tures in excess of 26 C are required to produce hurricane- 
force winds. Through these winds, the hurricane imparts 
stress to the ocean, which acts as a sink for the hurricane’s 
momentum. Interactions have been modeled in several ways. 

It is apparent that for an ocean-interacting hurricane model 
to be valid, the heat -momentum energy-exchange balance must 
be satisfied. 

Elsberry, Pearson, and Corgnati (1974) developed an axi- 
symmetric, ocean- interacting hurricane model which met the 
energy balance conditions. A description of the air-sea 
interaction within the boundary- layer of a hurricane was the 
output of and the primary purpose for the design of the model. 
For this purpose, a stationary, mature hurricane over a con- 
stant, isothermal ocean was assumed. The capability for in- 
teraction with a time-dependent ocean model was included in 
the design. 

One of the uses envisioned for the model described above 
was to investigate the cause of low sea-surface temperatures 
found in the wake of a hurricane. Based on limited observa- 
tions following tropical cyclones, wind mixing, convective 
mixing, and upwelling appear to be the main cooling mechan- 
isms. Jordan (1964) concluded that mixing was dominant. 
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Fischer (1958) mentioned upwelling as being important. 
Leipper (1967) found upwelling to be dominant near the eye 
with a transition to convective and/or mechanical mixing 
dominance at greater radii. The classical model by O'Brien 
and Reid (1967) and O'Brien (1967) simulated a stationary 
hurricane with constant forcing over a one-layer ocean. Re- 
sults of the model emphasized the importance of upwelling 
and downwelling. Gilbert (1974) modeled a stationary, time- 
dependent hurricane over a multi-layer ocean. Results of 
the Gilbert model indicate transitory inertial oscillations 
exert as much influence as the steady-state wind. 

Fraim (1973) developed a time -dependent , two-layer model 
for the ocean only based on the earlier mixed- layer models 
of Kraus and Turner (1967) and Denman (1973). The model was 
driven by the time -dependent , axi -symmetric , stationary 
hurricane model of Elsberry, et. at. (1974) . The axi-sym- 
metric assumption permitted representation of the hurricane 
and the ocean by a radial cross-section assumed to be the 
integrated average of all radial cross-sections. Fraim 
(1973) concluded that a combination of entrainment mixing 
and upwelling was the major factor in changing the mixed- 
layer temperature and depth. 

One of several logical extensions of the model was to 
use it to examine the effects on the ocean induced by a 
moving hurricane. One purpose for moving the storm was to 
consider Geisler's (1970) internal wake theory. The theory 
involves a comparison of the hurricane translational speed 
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(U) with the baroclinic wave speed (C) . For a slow moving 
storm (U/C < / 2 ) , Geisler postulated upwelling and resultant 
surface cooling. For a fast moving storm (U/C > 2) , Geisler 
postulated an internal wave wake at the interface between 
the mixed layer and thermocline. It should be emphasized 
that while Geisler's model and a similar model by O'Brien 
(1968) were primarily hydrodynamic, the model under study 
here was thermodynamic in nature. 

The present model was originally designed to be fully a 
three-dimensional model. However, due to computer time and 
storage constraints, the model was made to represent a ver- 
tical cross-section of the ocean along which the storm was 
moved. The storm was translated at several typical speeds 
(including U = 0) in order to determine the effects of move- 
ment on changes in the mixed- layer depth and temperature. 
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II. DESCRIPTION OF THE MODEL 



A. THE HURRICANE MODEL 

The hurricane model was a two-layer, time-dependent model 
by Elsberry, Pearson, and Corgnati (1974). By assuming the 
hurricane to be mature, an assumption of axi-symmetry could 
be made. A surface inflow layer and an upper outflow layer 
constituted the model. The tangential wind (Vg) profile at 
the top of the inflow layer was taken to be 

V 0 r 2 = constant (1) 

at radii outside the radius of maximum wind r^ . Inside r^ 
the winds were assumed to obey 

V 0 r * = constant. (2) 

Equations (1) and (2) were empirically derived from observa- 
tions by Riehl (1963) and supported by observations by Gray 
and Shea (1973). The maximum extent of the storm was de- 
fined to be the radius (r ) at which flow in the outflow 
layer turned from cyclonic to anti-cyclonic, i.e. Vq = 0. 

The radius of maximum wind was assumed to coincide with the 
eyewall radius, which was taken to be the radius at which 
the average of five adjacent values of equivalent potential 
temperature (6 ) was a maximum. Since the eye-wall has been 
observed to be a region of massive latent heat flux due to 
evaporation, the maximum 0 g was strongly tied to the sea- 
surfaCe temperature of the underlying water. An expression 



12 



for the maximum wind speed was determined by assuming a 
gradient balance between the wind field and the pressure 
field 
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p 8r r v e 
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(where f is the coriolis parameter) and using the relation 
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From r Q to 2t q the wind profile in (1) was modified to de- 
crease exponentially to 5 meters/second as a transition zone 
to realistic boundary conditions in order to allow downwelling. 

From observations by Riehl (1963) , the surface pressure 
at the eyewall was defined by the relation 



6p = 2.5660 (6) 

v 

where the 6 symbols indicate a deviation from p = 1005 milli- 
bars and 0 e =35O K. The radial pressure field was computed by 
integrating the gradient wind over the radius, using the wind 
speed from (1) and (2) . 

The atmospheric boundary-layer model consisted of a sur- 
face layer and a spiral layer above. Specification of the 
surface layer was by Monin-Obukhov similarity theory. Within 
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the surface layer, the surface roughness was a function of 
the stress at the surface. The spiral layer, a transition 
layer from radial to tangential flow, was specified by a 
modification of Ekman-layer theory. An iterative solution 
was used to match the two layers. Further details of the 
iterative solution can be found in Cardone (1-969) . 

A simple model for the hurricane was required to keep 
computer time and storage requirements to a minimum. Compar- 
atively simple specification of the boundary layer and lack 
of specification of upper layers enabled compliance with 
those requirements. The axi-symmetric assumption allowed 
representation of the storm by a radial cross-section. Again, 
computational considerations required such an assumption. A 
friction velocity for the computation of water velocities 
and kinetic energy gain, and the heat flux were provided to 
the ocean model by the hurricane model. A new radial distri- 
bution of mixed-layer temperature was predicted by the ocean 
model, and then provided to the hurricane model. This type 
of interaction can lead to storm intensification. However, 
for the purposes of this study, the hurricane maximum wind 
was bounded at 50 m/sec. This enabled isolation of changes 
in the ocean. 

B. THE OCEAN MODEL 

1 . Description of the Model 

The ocean model was a modification of that by Fraim 
(1973) , which was based on earlier mixed-layer models by 
Kraus and Turner (1967) and Denman (1972). The ocean was 
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assumed to have an initial mixed-layer depth of 30 m with a 
mixed-layer temperature of 30 C. The below-layer gradient 
was initially 0.1 C/meter to the top of an undisturbed layer 
defined by T = 20 C at Z = 130 meters. These values are 
typical of those found in the Gulf of Mexico and the Caribbean 
Sea in June, July, and August. The model was idealistic in 
that truly isothermal mixed layers and linear below-layer 
gradients are seldom found. Also, it would not be expected 
that the same thermal structure would exist over the whole 
ocean, or even any sizeable fraction of it. 

There were several important modifications of the- 
Fraim (1973) model. The Fraim model examined a radial cross- 
section of the ocean beneath a stationary hurricane. Equa- 
tions for the ocean model were written in R-6 co-ordinates. 

The present model examined a vertical cross-section of the 
ocean along which the diametrical cross-section of the hur- 
ricane was allowed to move. The diametrical cross-section 
consisted of the radial cross-section used by Fraim and the 
mirror image of that radial cross-section. This was made 
possible by the assumption of axi -symmetry in the hurricane 
model. The ocean equations were written in cartesian co- 
ordinates on a 200 X 3 grid (see Figure 1) with a 6 kilo- 
meter grid whereas Fraim used a 100 X 1 grid with a 3 
kilometer grid spacing. The time increment was 8 minutes 
compared to the 10 minute increment used by Fraim. A time 
increment of 8 minutes was used because it gave finer con- 
trol over available translational speeds. The movement of 
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BOUNDARY CONDITIONS: 

FOR OCEAN VALUES OF VELOCITY 
DEPTH. AND TEMPERATURE 
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Figure 1. Schematic of Hurricane Moving Along 200 X 3 Grid in the Open Ocean. 



the storm was achieved by moving the center of the storm to 
successive grid-points at time increments required to ap- 
proximate the desired speed. Another modification was the 
addition of a horizontal return flow within the thermocline 
layer extending down to a specified layer of no motion. This 
closed the ocean system except for the radial' advection, and 
thus improved the heat balance calculation for the stationary 
storm. Lastly, there was added to the model a convective 
adjustment for local convective instabilities. 

Several important assumptions were made about the 
manner in which the ocean and hurricane would interact. 

First, since the storm was assumed to be mature at initia- 
tion, the wind stress was imparted to the ocean exponential- 
ly so that after 9 hours the ocean felt 95 % of the stress. 
Second, some assumption had to be made as to what fraction 
of the stress would be used in inducing turbulence for en- 
trainment mixing. The balance of the stress would be used 
in the creation of horizontal currents. With no pertinent 
data available it was decided to let entrainment dominate 
initially. An exponential transition dependent on the mixed- 
layer depth shifted the dominance to currents as the mixed 
layer deepened. Third, it was assumed that ocean currents 
began in Ekman balance and adjusted immediately to changes 
in the stress pattern. This assumption eliminates inertial 
waves which must be resolved using smaller time steps. 

Fourth, in examining only a vertical cross-section it was 
assumed that there was no variation in temperature or ocean 
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velocity across the path of the hurricane. This simplifica- 
tion was made because the net. Ekman current should be normal 
to wind flowing inward at a specified angle, rather than 
being directed radially. Because of this, there was a re- 
sultant lack of symmetry across the path. This lack of sym- 
metry introduced instabilities into the model.. Fifth, it 
was assumed that the storm was wind dominated, implying that 
downward heat flux due to solar radiation was negligible. 

The ocean model was forced by kinetic energy input 
due to wind stress, upward sensible heat fluxes, and upward 
latent heat flux due to evaporation. Horizontal currents 
due to wind stress were computed. Resultant upwelling and 
downwelling due to divergent and convergent current regimes 
were then determined. A new mixed-layer depth was calculated 
through kinetic energy, heat flux, and advection considera- 
tions. Based on the new mixed-layer depth, a new mixed-layer 
temperature was calculated. The below-layer gradient and 
temperature at the top of the thermocline were then deter- 
mined. In this model the mixed- layer temperature and the 
temperature at the top of the thermocline can be different. 

If the temperature at the top of the thermocline was greater 
than the mixed-layer temperature a convective overturning 
process was initiated in order to stabilize the thermal 
structure. Following these calculations a radially averaged 
mixed-layer temperature was passed to the radial cross-section 
of the hurricane. 
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2 . Derivation of the Model Equations 

Following Denman (1973), the first law of thermo 
dynamics can be written 

-yz 



dT 

cTt 



YR*e 



P c 
K w p 



cnr~ 

K w p 



where Q^, = R*e the heat source term, 

= ocean water density, 
y = average extinction coefficient, 

0^ = specific heat at constant pressure, 

R* = solar radiation incident on the sea-surface, and 
z = depth in the ocean, positive downward. 

Expanding the total derivative in (7) yields 
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dz 



where w = is positive for downward motion. 
Time-averaging the turbulent form of (8) gives 
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3 z 



( 10 ) 



where the terms u'T' and v'T' were neglected because it was 
assumed that there was no systematic correlation between u' 
and v' and T' . 

Turning now to an examination of the boundary condi- 
tions, firstly at the ocean surface, z = 0, a heat loss at 
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the sea surface will induce convective overturning which 
cools the mixed layer. Therefore, the downward turbulent 



heat flux (w'T*) must equal the non- radiative heat loss. 



l . e . 


Qa 

( w ’ T ')o = • pV <“> 

K w p 


where 


= % * Qs- 


and Qg 


= upward latent heat flux due to evaporation, 


and Qg 


= upward sensible heat flux. 



At the base of the mixed layer downward turbulent heat flux 
is caused by wind-induced mixing and convective overturning. 
The turbulent heat flux must equal the heat transfer due to 
entrainment of water into the mixed layer. This balance is 
expressed by 
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where h 

0 


mixed-layer depth, 


T = 
s 


mixed-layer temperature, 


II 

H 


temperature at the top of the thermocline, and 


dh 

o _ 
dt 


the time variation of the mixed-layer due to wind 
mixing and convective overturning. 



For entrainment to occur (dh Q /dt) > 0, that is, there can be 
no shallowing of the mixed layer unless w is upward. Inte- 
grating (10) over the mixed layer and substituting the bound 
ary conditions gives 

v dT s Q t d h 0 Q a 

h o -at = p-?r h 0 • c - “ ) c T s • V - rc- (13) 

w p w p 
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Expanding the total derivatives, neglecting solar radiation, 
and solving for 3h Q /3t gives 
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(14) 



which is a predictive equation for deepening shallowing of 
the mixed layer in terms of the local change in mixed-layer 
temperature . 

To derive a similar predictive equation for the 
mixed-layer temperature Denman (1973) integrated (10) to an 
arbitrary depth z, and then over the mixed- layer depth 
yielding 
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where -/ 0 (w'T*)^ dz is the conversion of potential energy 
to kinetic energy by convection. If the turbulent kinetic 
energy is not varying in time, we may write 



W + G - D = 0 (16) 

where G = kinetic energy input from the wind, 

D = dissipation within the mixed layer, and 
h 

W = -j/* 0 (w 1 T ' ) dz 
o 

then substituting (16) into (15) , expanding the total deriva- 
tives , and solving for aT s /at yields 
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( 17 ) 
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which is the predictive equation for change in the mixed- 
layer temperature due to convection, entrainment mixing, or 
advection. Setting = 0 , and/or (G-D) = 0, and/or u, v, 
and w = 0 allows isolation of the cooling mechanisms for use 
individually and in combination. Substitution of (17) into 
(14) yields 
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(18) 



which is the predictive equation for changes in the mixed- 
layer depth as used in the present model. 

A partitioning procedure was used to divide the total 
available stress, t , from the boundary layer into one frac- 

3 . 

tion for the production of turbulence, x , and the remainder 
(t & - x Q ) , for the generation of wind-driven currents. 

h 

x = x EXP ( - pA- ) (19) 

o a v Cz J K ' 

o 

where x = total stress available from boundary- layer model, 

3 . 

h Q = mixed- layer depth, 
z q = roughness length, and 
C = a constant (10 1 *) (after Fraim 1973). 

Kraus and Turner (1967) assuming a value for x , 
specified the mechanical energy input (G) in equations (16), 
(17) , and (18) as 
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where g = gravitational acceleration, 

a = (1/p^) (dp^/dt) , coefficient of expansion 

p = air density, 

p = sea water density, 
w 

U* = atmospheric friction velocity, = /x /p , and 

3 3 

x = that part of atmospheric surface stress used to 
produce turbulent mixing. 

The dissipation (D) was set equal to zero. 

The calculation for the ocean velocities was by means 

of Ekman-spiral theory. The theory specifies a depth (Z ) 

at which the current direction has reversed from its surface 

direction, the current is EXP(-tt) of its surface value, and 

below which the wind has no direct influence. The velocity 

components of the ocean within the spiral layer are 

(x -x )cosa 

u= — - — EXP(-BZ) [cos (-BZ) -sin(-BZ) ] (21) 




(x -x )sina 

v= — - — ~ EXP(-BZ) [cos (-BZ) - sin (-BZ) ] (22) 

p (2f— ) ^ 
w p w 

where Z q = tt/B, specifies the Ekman depth, 

B = [f/2p/p J h , 

f = coriolis parameter, 

p = sea water density, 
w 

p = vertical eddy coefficient of viscosity, 
z = depth, and 

a = atmospheric inflow angle. 
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Integration of (21) and (22) over the Ekman depth (Z ) gives 

V 

the net mass transport in the Ekman layer 



M x = (x a - T 0 ) cosa/ f (23) 

M y = ^ T a " T 0 ) sina / £ * (24) 



From these mass transport equations it is possible to calculate 
the mean ocean velocities in the mixed layer 
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The horizontal velocities in the layer are related to the ver- 
tical velocity w by the continuity equation 
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Equations (17) and (18) give two equations for three 

unknowns (T ,h ,T, ). A return flow is calculated for the 
v s’ o’ rr 

layer between 130 meters and the new mixed- layer depth. It 
is assumed to have the same mass transport as the opposing 
current in the mixed layer, and to be centered on a point 
2/3 of the distance from 130 meters to h . The below-layer 
gradient is adjusted by the return flow, upwelling, and down- 
welling. The new T^ is the temperature at the intersection 
of the new below-layer gradient and the new mixed- layer 
depth. This relation is given by 




3T 

3z 



) (Z 



B 



h ) 
o' 



24 



where 3T/3z = is the below-layer gradient, 



T r = temperature at top of the undisturbed layer 
= 20 C, and 

Zg = depth of top of undistirbued layer = 130 m. 



25 



III. DISCUSSION OF RESULTS 



It was the intent of this study to use a sea-air inter- 
acting numerical model to simulate the effects upon the 
ocean of a hurricane moving across it. The reasons for con- 
ducting the study were: 1) to investigate the cause of low 

sea-surface temperatures in the wake of a hurricane, and 2) 
to determine the effects on the ocean for different speeds 
of storm movement. 

Changes to the mixed-layer depth are controlled by Equa- 
tions (17) and (18). Convective mixing is represented by 
the first term in both equations. Entrainment mixing in- 
volves the second term in both equations. Advection is 
described by the remaining terms in the equations. Accord- 
ing to the two equations entrainment mixing and convection 
can only cool and deepen the mixed layer. It has been as- 
sumed that is positive by assuming a wind-dominated 
regime. On the other hand, since ocean velocity values can 
be both negative and positive, the advection terms can both 
warm and cool, shallow and deepen the mixed layer. 

This discussion is divided into three parts: 1) a com- 

parison with the results of Fraim (1973) , 2) a discussion of 
the relative importance in this model of the mechanisms dis- 
cussed above, and 3) a discussion of the effect of the hurri 
cane translational speed on changes in the ocean. 
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A. COMPARISON WITH FRAIM'S RESULTS 



Before attempting to move the hurricane, it was necessary 
to establish that the stationary version of the present model 
yielded reasonable results. For the initial experiment, the 
200 x 3 grid (Figure 1) was defined with a 3 km grid spacing 
as was the 100 x 1 grid in Fraim's model. Profiles of values 
at 12 hours for the present model are shown in Figure 2. A 
comparison of maximum and minimum values for Fraim's model 
and the present model are listed in Table I. 

The first point of comparison between the two studies is 

the hurricane portion of the model. The initial conditions 

for both hurricanes were the same: air temperature T a = 29C, 

radius of maximum wind r^ = 27km, and maximum tangential wind 

V_ = 30m/sec. There were some differences. Fraim did not 
u . 

1 

constrain the maximum tangential wind or the radius of maxi- 
mum wind. In the present model, for this experiment, the 
maximum tangential wind was forced to increase sinsuoidally 
from 30m/sec to 50m/sec at 12 hours. Also, a limit of 
r^ < 60 km was set on the radius of maximum wind. 

The tangential wind profile (bottom. Figure 2) is char- 
acterized by a rapid increase with radius to its maximum 
value at r Q = 60 km (Equation (2)), a slow decline to its 
value at r Q (Equation (1)), and a more rapid decline to the 
outer boundary value. From r Q to 2r Q the wind is allowed to 
drop off to 5m/sec. At r Q the change in the wind profile is 
abrupt and therefore easily detected. The radius r Q tends 
to be the limit of applicability of the assumption of axi- 
symmgtry at upper levels, although the assumption is good to 
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o 




TIME 12-0 HOURS 



Figure 2. 12-Hour Predicted Values of Tangential Wind, 

Ocean Radial Velocity, and Ocean Vertical Velocity 
for the Stationary Diametrical Hurricane Model on” 
a 3 km Grid. 
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Table I. Comparison Between Some 12-Hour Values of the 

Fraim Radial Stationary Hurricane Model and the 
Diametrical Stationary Hurricane Model. 



Max ORV 
Max DW 
Max UW 
Min MLD 
Max MLD 
Min MLT 
Max MLT 
Max TW 
RMW 

Area average ocean 
heat loss 

Min TT 

Max TT 



Fraim Model 

.48 

.03 

3.35 

4.31 

51.23 

25.26 

30.43 

55.58 

42 

930 

19.80 

29.99 



ORV - ocean radial velocity (m/sec) 
DW, UW - downwelling, upwelling (m/hr) 
MLD - mixed-layer depth (m) 

MLT - mixed-layer temperature (C) 

TW - tangential wind (m/sec) 

RMW - radius of maximum wind (km) 

TT - thermocline temperature (C) 



Diametrical Model 

.28 

.21 

1.88 

21.23 

58.20 

28.84 

30.00 

50 

60 

371 

26.90 

30 



Heat loss (cal/cm 2 ) 
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a greater radius at lower levels. Because r Q expands with 

intensification (r = 225 km at 12 hours vs. 150 km at 0 

v o 

hours), the wind velocity value at the outer boundary of 
the hurricane model has increased from 5m/sec to 12.5m/sec. 
The maximum tangential wind (Table I) was about 101 greater 
in Fraim's model than that of the present model due to lack 
of constraint on the maximum wind value. 

The ocean radial velocity (actually the x-component of 
the horizontal velocity) (middle , Figure 2) responds directly 
to the stress imparted by the tangential wind (Equations 
(21) and (22)). Ocean velocity values are calculated by 
integrating the x-component of the Ekman current over the 
mixed-layer depth. Fraim integrated the mass transport over 
a shallower depth thus giving a larger velocity (Fraim's 
results , Table I) . 

As the ocean radial velocity directly follows the tan- 
gential wind, the ocean vertical velocity (top, Figure 2) 
follows the gradient of the radial velocity. The result is 
that the maximum upwelling occurs at the point of maximum 
increase of ocean radial velocity, i.e., between the center 
of the storm and the radius of. maximum wind. Maximum down- 
welling occurs at r Q where the ocean radial velocity drops 
sharply. Fraim's model upwelling (Table I) is stronger be- 
cause the radius of maximum wind was smaller. The Fraim 
model had no major downwelling because tangential wind pro- 
file was not altered in the r to 2r region. 

The radius of the maximum mixed layer depth, minimum 
mixed-layer temperature, and minimum thermocline temperature, 
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(Figure 3) is the same, and coincides with the radius of 
maximum wind. Coinciding with the radius of maximum upwell- 
ing are the minimum mixed layer depth, and pips of warming 
in the mixed-layer (T g ) and thermocline (T^) temperatures. 

The pips appear to be due to a lack of cooling due to upwell- 
ing. The minimum mixed-layer temperature (Table I) is much 
smaller in Fraim's model because it did not have the thermo- 
cline layer return flow discussed in Chapter II B. 1. Lack 
of a return flow of water allowed strong cooling of the 
thermocline layer in the Fraim model which created a large 
T g -T^ difference (Table I) and retarded deepening of the 
mixed layer. Consequently, the area- averaged heat loss in 
Fraim's model was more than twice the loss in the present 
model (see Table I). Most of the heat loss is concentrated 
within 100 km of the center. Advection of water both hori- 
zontally and vertically acts to cool or warm the thermo- 
cline layer. Warming in the upper portion of this layer 
increases the below-layer gradient while cooling decreases 
the gradient. A change in the below-layer gradient is thus 
associated with the temperature at the top of the thermo- 
cline, that is, the thermocline temperature, since the tem- 
perature at the level of no-motion is held constant. 

Stronger upwelling in Fraim's model (Table I) accounts for 
the shallower mixed-layer depth minimum and seems to indicate a 
higher maximum mixed-layer temperature. The comparison be- 
tween the models shows the results are similar, however, due 
to the differences mentioned above the present model yielded 
more, reasonable results. 
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MIXED LR V ER MIXED LRVER THERMOCLINE 

DEPTH (M) TEMP. (C) TEMP. .fC) 




Figure 3. 12-Hour Predicted Values of Mixed-Layer Depth, 

Mixed-Layer Temperature, and Thermocline Tempera- 
ture for the Stationary Diametrical Hurricane 
Model on a 3 km Grid. 
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B. DISCUSSION OF COOLING MECHANISMS BENEATH A MOVING STORM 



As the comparison of the present model with Fraim's 
was considered satisfactory, it was decided that the results 
warranted progressing to a moving storm. The storm was al- 
lowed to interact with the ocean until a maximum wind of 
50 m/sec was achieved so as to allow examination of the 
ocean's response to the same storm. A new limit of r^ < 45 
km was placed on the radius of maximum wind. The 200 x 3 
grid spacing was expanded to 6 km giving a 12 km by 1200 km 
rectangle in the open ocean. Starting point for the center 
of the storm was the 900 km point (Figure 8) . Movement of 
the storm was from right to left at 1 grid increment per 2, 

3, 5, and 7 time increments, or 6 km per 16, 24, 40, and 56 
minutes which corresponds to 11.84, 7.89, 4.75, and 3.39 
knots respectively. 

A time cross-section of the 4.75 knot storm was drawn 
for the mixed-layer temperature response (Figure 4, Table II) 
and mixed-layer depth response (Figure 5, Table III) of the 
ocean to the mechanisms indicated in the figures. The 750 
km point was chosen for this section because during the 30 
hour forecast period it experienced passage of the most in- 
teresting parts of the storm i.e. the radii of maximum wind, 
the radii of maximum upwelling, and the no-stress zone of 
the eye. Figure 8 can be used to see the relative locations 
of these zones. 

As shown in Figure 4, advection alone has virtually no 
impact on the mixed-layer temperature, though in the absence 
of deepening (Figure 5) shows a slight warming as mentioned 
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Table II. Actual Values of Mixed-Layer Temperature for Key Points of Passage at the 750 km 
Point of the 4.75 Knot Storm. 
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Table I for abbreviations 



in Chap. III. A. Entrainment mixing and convection appear 
to dominate outside the radii of maximum upwelling with the 
maximum effect being inside the radii of maximum wind. Of 
the two, entrainment mixing is the stronger cooling mech- 
anism. It is interesting to note that convection alone con- 
tributes to the small amount of cooling within the radii of 
maximum upwelling. All mechanisms combined appear to con- 
tribute to cooling between the radii of maximum upwelling, 
although no changes occur during passage of the eye. It ap- 
pears that horizontal ocean velocities have little direct 
effect on the mixed layer, but that vertical velocities have 
a significant effect. 

As shown in Figure 5, advection alone has a small effect 
in changing the mixed layer. Again, entrainment mixing and 
convection appear to dominate outside the radii of maximum 
upwelling with the maximum effect being inside the radii of 
maximum wind. Of the two, entrainment is definitely dominant 
in deepening the mixed layer, although inside the radii of 
maximum upwelling convection is the only mechanism that 
deepens the mixed layer. Inside the radii of maximum upwell- 
ing, advection alone is responsible for shallowing the mixed 
layer. In conjunction with entrainment mixing and convection, 
the effect of advection increases. 

One cannot completely separate the mechanisms since there 
is a complex interaction between them in the ocean. However, 
Leipper (1967) interpreted baththermograph data before and 
after passage of the hurricane Hilda in the Gulf of Mexico 
during 1964 as follows: 1) cold water was upwelled due to 
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divergence in the region near the center of the storm, and 
2) warm water farthest from the storm path gave evidence of 
having been mixed by surface cooling or by mechanical mixing 
or both. As shov/n in Figure 9, the lowest mixed-layer tem- 
peratures are in the vicinity of the storm center. The 
warmer water toward the boundary of the storm' (Figure 8) 
shows evidence of mixing (implied by mixed layer deepening) . 

A set of profiles of temperature change from the initial 
temperatures in the ocean due to hurricane passage is shown 
in Figure 6. The profiles show cooling down to 30 m (except 
the profile for advection only) then warming down to the 
mixed-layer depth. Entrainment of colder thermocline water 
cools the mixed layer and warms the thermocline layer to 
form a new mixed layer. At the mixed-layer depth in the 
model there is a discontinuity denoted earlier by T g - T^. 
Below the mixed layer only those mechanisms involving advec- 
tion show distinct cooling. Return flow in the thermocline 
layer and upwelling are responsible for the cooling. The 
750 km point is 120 km behind the storm center at 30 hours. 
As a result Figure 6 agrees with Federov's (1973) Type A 
thermal structure change profile of two layers of cooling 
separated by a layer of warming. 

The effect of the return flow in the thermocline layer 
can also be seen in Figure 7. Entrainment mixing and con- 
vection singly and together have an insignificant effect on 
the below-layer gradient. The rapid decrease of the below- 
layer gradient for advection alone shows that for intense 
upwelling there is strong cooling in the upper part of the 
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Figure 6. Change in Thermal Structure at the 750 km Point 
After Passage of the 4.75 Knot Storm. (The 
initial conditions were MLD = 30 m, MLT = 30 C 
with a linear temperature decrease to a value 
of 20 C at 130 m.) 
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thermocline layer indicating a more intense upwelling than 
for advection alone. 

C. EFFECT OF STORM MOVEMENT ON OCEAN RESPONSE 

The second purpose listed for the study was to determine 
the effect of movement of the storm on the ocean's response 
to the storm. In Table IV is listed some maximum and minimum 
values of key parameters for several storm translational 
speeds . 

Changes in the mixed layer due to storm movement are 
comparatively easy to explain. Slower-moving storms have 
more time at a given location for entrainment mixing and con- 
vection to cool and deepen the mixed layer. For the sta- 
tionary storm the downwelling maximum nearest the storm 
center (see Figure 10) has added to the deepening but per- 
haps retarded the cooling. The anomalous spikes of down- 
welling (for instance at the 900 km point in Figure 10) are 
due to the sharp drop in the tangential wind at the storm 
boundary and the mixed-layer depth spikes generated due to 
the stationarity of the storm. 

Increasing heat flux to the storm with increasing storm 
speed (Table IV) is due to fast moving storms extracting more 
heat because of continuously encountering uncooled 30 C 
water. Radial heat flux (Table IV) is due to heat flow in 
currents out through the mixed layer and in through the 
thermocline layer. For a stationary storm advection supplies 
more heat than it removes. This seems to be the source of 
warm water for the mixed-layer temperature pips (Figure 11). 
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TRNGENTIRL OCEPN RRDIPL 0CERN UERTICRL 

WIND (M/SEC) UEL. CM/SEC) UEL- (IVH0UR) 

Z5 50 -0.1 0-0 0.1 -1.5 0-?5 



o 




TIME -- 18- D HOURS 

Figure 8. 18-Hour Predicted Values of Tangential Wind, Ocean 
Radial Velocity, and Ocean Vertical Velocity for a 
Hurricane Moving at 4.75 Knots. 
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MIXED LAYER MIXED LAYER THERM0CLINE 

DEPTH (M! TEMP. (Cl TEMP. (C) 



on 





TIME -- 18-0 HOURS 

Figure 9. 18-Hour Predicted Values of Mixed-Layer Depth, 

Mixed-Layer Temperature, and Thermocline Tempera- 
ture for a Hurricane Moving at 4.75 Knots. 
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Figure 10. Same as Figure 8 Except for a Stationary Hurricane. 
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TIME -- 18. G HOURS 

Figure 11. Same as Figure 9 Except for a Stationary Hurricane. 
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Table IV. 



Comparison of Parameters for Several Storm 
Translational Speeds at t = 18 Hours. 





Stat 


3.39 
kts . 


4.75 
kts . 


7.89 
kts . 


11.84 
kts . 


max ORV 


.27 


.28 


.29 


.29 


.29 


max DW 


.48 


.32 


.28 


.27 ' 


.21 


max UW 


1.18 


3.10 


2.98 


2.78 


2.60 


min MLD 


17.18 


30 


30 


30 


30 


max MLD 


72.60 


61.62 


57.91 


54.78 


53.94 


min MLT 


28.85 


28.86 


28.88 


29.06 


29.25 


max MLT 


30.09 


30 


30 


30 


30 


area average 
ocean heat 
loss 


- 29 


265 


325 


352 


335 


radial heat 
flux 


-146 


58 


89 


103 


92 


heat lost to 
storm 


156 


161 


163 


168 


174 



*see Table I for abbreviations 
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The apparent reversal of radial heat flux values for the two 
fastest storms in unexplained. 

A sustained upwelling for 18 hours duration would account 
for the shallow minimum mixed-layer depth (Table IV) . By the 
equation for continuity (Equation (27)) the integrated ver- 
tical velocity is equal to the horizontal velocity divergence 
times the integrated depth. Therefore the maximum upwelling 
value for the stationary storm is small because of the mixed- 
layer depth in which it occurs (Figures 10 and 11) . For a 
moving storm the leading upwelling zone the upwelling is more 
intense because it occurs in the deepest mixed layer (Figures 
8 and 9) . Maximum downwelling intensities are governed by 
mixed-layer depth also (Figures 8, 9, 10 and 11). Maximum 
values for the ocean radial velocity are much the same be- 
cause the mean mixed-layer depths at the radius of maximum 
ocean radial velocity are nearly the same. The ocean radial 
velocity profiles are altered slightly due to mixed-layer 
depth considerations (Figures 8 and 10) . 
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IV. SUMMARY AND CONCLUSION 



It was the intent of this study to use a numerical model 
to simulate the effects of a hurricane moving across the 
ocean. In such numerical models the inherent limitations in 
elude the finite differencing scheme as well as the initial 
and boundary conditions. Imposed limitations involve assump 
tions associated with applying the model. 

For this study a simple model was used for both the hur- 
ricane and the ocean. A bound was set on the hurricane 
maximum wind in order to standardize the storm input to the 
ocean. There were several key assumptions about the way in 
which the ocean and storm interact: 1) wind stress was par- 

titioned into making currents and turbulence based on the 
mixed-layer depth, 2) ocean currents were assumed to be in 
Ekman balance and to respond immediately to changes in the 
stress pattern, 3) it was assumed that there were no gradi- 
ents of temperature or ocean velocity across the path of the 
storm, 4) the storm was assumed to be wind-dominated, and 
5) since the storm was mature at initiation the total stress 
imparted exponentially until 95% of the total stress was 
imparted to the ocean at 9 hours. 

In answer to the two main questions of interest: 1) 

cooling of the ocean appears to be primarily due to entrain- 
ment and convection except in intense upwelling within the 
region of maximum wind. In this region all mechanisms con- 
tribute to cooling of the mixed layer, 2) slow moving storms 
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cool and deepen the mixed layer more than do fast moving 
storms. Waves in the mixed-layer depth due to internal 
waves on the mixed-layer thermocline layer interface postu- 
lated by Geisler (1970) were not seen. This is due to the 
thermodynamic nature of the model under study here compared 
to the hydrodynamic nature of Geisler’s model. 

In summary, the results are satisfactory within the limits 
of the model. The logical course to follow is to reduce the 
number of assumptions used. In a model with as many simul- 
taneous interactions as this one, the task is one to be taken 
slowly, in small increments. 
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